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Schmid  and  Frye  (1976)  presented  a  plan  for  rating  the 
susceptibility  of  Engelmann  spruce-subalpine  fir  {Picea 
engelmannii  Parry  and  Abies  lasiocarpa  (Nutt.)  Hooker) 
Stands  to  spruce  beetle  (Dendroctonus  rufipennis  (Kirby)) 
infestation.  The  plan  ranks  the  susceptibility  of  the  stand 
based  on  its  physiographic  location,  basal  area,  average 
diameter  of  the  spruce,  and  percentage  of  spruce  in  the 
overstory. 

During  the  past  5  years,  inventory  projects  in  the  Rocky 
Mountain  Region  of  the  USDA  Forest  Service  have  been 
gathering  intensive  data  on  spruce-fir  stands  by  the 
methods  outlined  in  the  1980  edition  of  the  Region's 
"Silvicultural  Examination  and  Prescription  Handbook" 
(FSH  2409.260).  Substantieil  information  is  gathered  for 
each  stand,  some  of  which  is  used  in  the  Schmid-Frye 
system.  The  data  are  recorded  on  field  forms  and  eventually 
are  transferred  to  magnetic  tapes  which  are  stored  in  the 
project's  headquarters. 

Since  its  origin,  the  rating  system  has  received  limited  use 
in  the  southern  and  central  Rockies.  To  increase  its 
usefulness  to  forest  managers,  a  computer  program  was 
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written  to  calculate  stand  susceptibility  from  the  inventory 
data.  The  program  reads  the  inventory  file  for  the  stand, 
extracts  the  essential  information,  computes  the  stand 
rating,  and  prints  it  out,  identifying  the  stand  by  the  location 
data  commonly  recorded. 

This  program  may  give  a  slightly  different  rating  for  each 
stand  from  stage  II  inventory  data  than  it  would  if  the  data 
were  gathered  in  the  same  manner  as  was  used  to  develop 
the  rating  system.  In  Schmid  and  Frye  (1976),  the  limits  for 
each  category  in  the  basal  area,  diameter,  and  percent 
spruce  characteristics  were  based  on  data  derived  from 
fixed  radius  plots.  Stage  II  data  are  derived  from  variable 
radius  plots  and  different  basal  area  factors.  How  this  will 
affect  the  stand  rating  is  not  known  at  this  time,  but  no 
significant  change  is  assumed. 

The  program  is  compatible  with  USDA  Forest  Service, 
Rocky  Mountain  Region  program  processing  data  collected 
from  the  stage  11  inventory.  The  program  could  be 
integrated  into  the  master  program  of  the  Region  to  run  as  a 
subroutine  so  that  the  stand  rating  for  spruce  beetles  would 
be  automatically  included  in  the  printout  with  the  other 
stand  information  regularly  presented.  There  is  also  the 
possibility  of  a  computer-produced,  differentially  shaded  or 
colored,  overlay  map  representing  a  mosaic  of  standing 
ratings  for  the  different  stands  in  the  spruce-fir  forest.  If  the 
Schmid-Frye  rating  system  is  adopted  with  modifications 
for  regioneJ  conditions  in  other  western  regions  of  the 
USDA  Forest  Service,  this  program  can  also  be  easily 
modified  for  use  with  it. 
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Appendix  I 


Once  data  for  a  specific  stand  has  been  isolated  from  the 
data  file  for  all  the  stands  using  an  add  command  in  the 
Univac  1100  text  editing  routine  from  the  Univac  Manual 
UP-8723,  program  processing  requires  no  additional  input. 
The  inventory  file  for  the  stand  is  searched  by  the  program 
for  information  required  for  computation  of  stand  spruce 
beetle  susceptibility  index.  The  program  calls  USDA  Forest 
Service  subroutine  ESSITE  for  computation  of  site  index 
based  on  Alexander  (1967).  On  output,  the  program  lists: 

1.  Alexander's  site  index 

2.  Total  number  of  trees  processed 

3.  Number  of  sample  trees  with  d.b.h.  greater  than  2.6 

inches 


4.  Number  of  sample  spruce  with  d.b.h.  greater  than 

10.0  inches 

5.  Stand  risk  rating  from  Schmid  and  Frye  (1976) 


Schmid /Frye  Value 

4-5 
6 

7-9 
10 
11-12 

See  example  problem 


Program  Output 

Low 
Low  Medium 

Medium 
Medium  High 
High 


Appendix  II 

Program  listing 


C   PROGRAM   TO  ESTIMATE   SUSCEPTIBILITY  OF   A  SPRUCE  STAND  TO 

C   SPRUCE  BEETLE  ATTACK 

C 

C   DATA  INPUT   IS  STANOARO   REGION-2  STAGE   I/II   STAND  DATA 

C   IT   IS  ASSUMED  THAT  DATA   IS  FOR   SPRUCE  STANDS  ONLY 

C   AS  MANY  SPRUCE   STANDS  AS  NECESSARY  CAN  BE   IMPUTED  AT   ONE  TIME 

C**«*»*  USE   UNIT    7  FOR   DATA  INPUT 


DIMENSION  X(4>.IASCI(2t3) » AS  I ( 2 0 1 1 1 TI T L ( 12 > 

DATA    IASCI/4HHIGHtlH   f4HME3   , 1HH I6H » 4HMEDI ♦ 2HUM t^HLOU  ,3HME0» 
13HL0UflH  / 

C  ... 

C  ...   TDRH  =   TOTAL   LIVE  OBH   .ST.   2.6  IN 

C  ...   TSDRH  ~   TOTAL   LIVE   SPRUCE   .GT.   2.6  IN 

C  ...   TC  =   TOTAL   OBH   IN  CANOPY   .5T.    2.6  IN 

C  ...   TSC  =   TOTAL  SPRUCE   IN  CANOPY   .GT.   2.6  IN 

C  ...   SD8H10   =   TOTAL   SPRUCE   .GT.    IN  REGARDLESS  OF  POSITION 

C  ...   XBAI   IS  BASIAL   AREA  INDEX 

C  ...   XHI   =   TREE  HEIGHT 

C  ...  NUM=NUHBER  OF   ADDITIONAL  STEMS 

C  ...   XAGE   =   TREE  AGE 

C  ... 

C  ... 

C  ...   ICNT  =  COUNTER  FOR  TOTAL   TREES  PROCESSED 

C  •••   ICNT2   =  COUNTER  FFR   TREES   OVER   2.6   IN.  DBH 

C  ...   ICNT3   =  COUNTER  FFR  SPRUCE  OVER  10.  DBH 

C  ...   NPLTS  =   COUNTER  FOR  NUMBER  OF  PLOTS 

C  •••   IPLTSl  =   SECOND  PLOT  COUNTER 
C  ISITC  =  COUNTER  FOR  SITE  TREES 

C  ... 
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C   ...   RSK  =   SITE  INDEX 
TDBH=0- 
TSD8H=0. 
TC=0. 
TSC=0. 
SDBH10=0  . 
IUNIT=7 

C 

C  READ   HEADER  CARD 

197  CONTINUE 
NPLTS=0 
IPLTS1=0 
ICNT=0 
ICNT2=0 
ICNT3=0 
ISITC=0 
IEND=0 

READ(IUNITil98f EN0  =  98)   IDXl , IDX2 » CI  TIT L C I ) f I  =  1 1 1 2 » 

198  F0RMAT{2I6,lX»llA6tAl) 
IUNIT=7 

C   ...   DUMMY   READ  LOOP  TO  GET  AT  DATA 

20  READ(IUNIT»199»END=g8)   I CT • I  ASCI  2 f I ASC II 
C  URITE(S,ig9)ICT,IASCI2f lASCIl 

195     F0RMAT(12X»I1 tl6X»A2taXfAl> 
IF(ICT.NE.2)G0  TO  21 
DECODE < 1  ,50  It lASCIllIRSK 

501  FCRMATCIl) 
C   •  •  • 

C   III   XBAI   IS  BASAL   AREA   INDEX  FROH  INPUT 

555     FORMATI10X,"IRSK,XaAI",I5f 10X,IPE14.7I 
C    •  •  • 

DECODE (2 ,502, I ASCI2) XBAI 
C  WRTTE(6,555)IRSK,XBAI 

502  F0RMAT(F2.0) 
GO  TO  20 

21  CONTINUE 

IF(ICT.NE.8)G0   TO  20 
C    •••   READ   DATA  CARDS 

1  READf lUNITf 101 »END=98tERR=99>ISITl,ISIT2tICTf IPLTS«TTH»ITSP»OBHf 
IXHIffNUKf XAGEflTPOS 
101  FORMAT(2I6tTl,I3t7X»I2tI3tF3.1tF3.0»I2»5XfF3.0»10X»Il> 

C   CHECK  FOR  NEW  STANO 

IFCICT.EQ.l)   GO  TO  505 
C  ISIT  =   SITE  ID 

C   •••   ICT=  CARD  TYPE 

C  ITH  =  TREE  HISTORY  —  "Ol"  IS  LIVE  TREE  —   "gS*   IS  SITE  TREE 

C  ITSP  =   TREE  SPECIES 

C   •••   DBH  =  DBH 

C    •••   IRSK  =  PHYSIOLOGICAL  LOCATION  ■IRSK=8   .THEN.   HIGH  RISK 

C   III   BRANCH  TO  READ   IF  NOT  LIVE  TREE 

IF(IPLTS.NE.IPLTS1)NPLTS=NPLTS*1 

IPLTS1=IFLTS 
C   •••  CHECK  FOR  SITE  TREE 

IFCITH.EG.DGO  TO  2 

IF(ITH.NE.99)G0  TO  1 

C   •  •  • 

C   •••   TREE    IS  SITE  TREE 
C   ...   ASSUME   TREE  IS  SPRUCE 
C   ...  FROM  FOREST  SERVICE  SUBROUTINE 
ISITC=ISITC*1 

CALL  ESSITE(XHI,XAGE»ASI<ISITC) ) 
60  TO  1 
2  CONTINUE 

NUM=NUM  -  1 
IF(NUM.LT.O)  NUM=0 
ICNT=ICNT  ♦  1   ♦  NUM 
NUM=0 

C   ...   DISREGARD  TREE   IF  LESS  THAN  2.6  IN. 

IF(DBH.LT.2.6)G0  TO  1 

ICNT2=ICNT2*1 

DDDH=(DBH/2.)/12. 

DDH=DDDH*DDDH*3.2425962 
C  yRITE(6,500 )OBH,DDOH 

5C0  F0RMAT(10X,2(1PE14.7>> 

TDBH=T08H*DD0H 

IF(ITSP.GT.90 )TSC=TSC*ODOH 
6  IF(DBH.LT.10. )G0  TO  1 
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IF(ITSP.LT.90)SO  TO  1 
ICNT3=TCNT3*1 
SOBHie=SDBH10*08H 
60  TO  1 
38  IEND=1 

60  TO  505 
^=  CONTINUE 

C   CHECK  FOR  ERROR   IN  OATA  FORMAT 

IF(ICT.EQ.l)   GO  TO  505 
WRITEf6,504)  ISIT1»ISIT2 
504  FORHAK"   THERE  HAS  SEEN  AN  ERROR   IN  READING  THE  OATA  FOR" 
1"   STAND     "faief/f"       RUN  TERMINATED") 
STOP 
5C5  CONTINUE 

IF(IPLTS1«EQ.NPLTS)G0  TO  7 
C  URITE(6»503)NPLTS 
?C3     F0RHAT(10Xt"£RROR  —  PLOT  NUMBER  AND  COUNT  NOT  EQUAL"* 

1/flOXf "NIMBER  OF  PLOTS  COMPUTED  AS  =  "fl^l 
?  CONTINUE 

C  ...  COMPUTE  AVERAGE  DIAMETER  OF  LIVE  SPRUCE  ABOVE  10.  IN.  DBH 
C   ...  —  ALL  LIVE   SPRUCE  TREES  USED*   BOTH  CANOPY  AND  UNDERSTCRY 

IF(ICNT3.GT.0>GO  TO  10 

WRITEC6»300) 

2C0     FORMAKIOX, "STAND  CONTAINS  ZERO   SPRUCE   ABOVE   10.   IN.  DBH"»/ 
If 10Xt"VARIA8LE  XC2I   SET  TO  ZERO"! 
X(2>=0. 
60  TO  19 
10  X1=FL0AT(ICNT3) 
C   ...   *  *  NOTE  XBAI   IS  BASAL  AREA   INDEX  •  • 
X(21=SDBH10/X1 
19  CONTINUE 

C  ...  COMPUTE  TOTAL  STAND  DBH  REGAROLFISS  OF  SPECIES  OR  TREE  POSITION 
C   ...   TREES   MUST  BE  GREATER  THAN   2.6  DBH   TO  BE  INCLUDED 

X1=ICNT2 

X2=NPLTS 

X<3)=€X1/X2)*XBAI 
C   ...   COMPUTE   PERCENT  SPRUCE  IN  CANOPY  -  TREES  .GT.   2.6  DBH  ONLY 

X<4)=100.*(TSC/TOBH1 
C   ...   CONVERT   RAy  VALUES  TO  RISK  VALUES 
C  URITE(6*S00)XCl)fXC2)«Xf 3)tX(4) 

=  00  FORMAT(10X»4(1PE14.7M 
C   ...   SITE  INDEX 

IFCISITC.LE.OJGO  TO  600 
C  URITECS,*)"ISITC  =  "tlSITC 

X1=ISITC 

RSK=0. 

00  40  I=1,ISITC 
C  URITE(6,*)-ASI ("♦It")  =  "tASICD 

RSK=RSK-»ASICII 
40  CONTINUE 

RSK=RSK/X1 
X<11=2. 

IFfRSK.LT.80.)X(ll=l. 
C   ...  CHECK  FOR  CREEK  BOTTOM  OR  DRAW 

IF(IRSK.EQ.8>X(1)=3. 
C   ...   BEGIN  OUTPUT 

URITEC6,603)( ITITL(I) tl=l*12) 
e03     F0RMAT(//,10X»"SPRUCE  BEETLE  STAND  RATING  FOR  "vllASvAll 

URITE(6t601)RSK 
€01     F0RMATC10X»"ALEXANDERS  SITE   INDEX  =  •tF10.4> 

X1=X(2) 

X(2)=3. 

IF(X1.LT.16.>X(2)=2. 

IFIX1.LT.12.)XC2>=1. 

X1=X<3> 

X«3>=3. 

IF(X1.LT.150. )X(3>=2. 
IFCXl.LT .100. )X(3)=1. 
X1=X(4) 
X(4)=3. 

IF(X1.LT.65.)XC4)=2. 
IF<X1.LT.50. )X(4)=1. 
C   ...  COMPUTE   RISK  SCORE 
XRSK-0. 
DO  30  I=l»4 
XRSK=XRSK*X(I) 
30  CONTINUE 
I0UT=1 
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IFCXRSK.LT.ll.  )I0UT-2 
IF(SRSK.LT. 10. ) I0UT=3 
IFCXRSK.lt. 7. )IOUT=% 
IF(XRSK.LT.6. ) I0UT=5 


201 


IFCXRSK.lt. 7. )IOUT=% 
IF(XRSK.LT.6. ) I0UT=5 

WRITE(6f 301 >ICNTtICNT2tICNT3»CIASClCI»I0UT>  fl  =  l t2> 
FORMATf IGXt//, 
110X»"TOTAL  TREES  PROCESSED  IN  STAND  =  •»I3f 
2/10X»*TREES  PROCESSED  WITH  DBH  X2.6   IN.  =  "tlSt/ 
310X»"TREES  PROCESSED  WITH  OBHXIO.O   IN.   -  "tlS*/// 
%10X,«STANO  RISK  RATING  FOR   SPRUCE  BEETLE  IS  ■»2A4», 
 CHECK  TO   SEE   IF  END  OF  DATA 


■»2A4»//I 

 CHECK  TO   SEE   IF  END  OF  DATA 

IFIIEND.NE.O)  STOP 

 REREAD   TITLE  CARD  FOR  NEW   STAND  AND  BEGIN  PROCESSING   N«:U  STAND 

TiimTTz.^n 


IUNIT=30 
GO  TO  197 


600  CONTINUE 

WRITE(69602) 

602     F0RHAT(10X»"NO  SITE  TREES  —   PROGRAM  ABORT") 
STOP 
END 


C 
C 
C 
C 


12=1 

IFCA6CT)-TEH)  10f20,20 
10  CONTINUE 
20  00  30  J=2«ll 

J1=J-1 

J2=J 

HROWC J)  =  CHTSI(I1»J)»CTEN-AG< 12 » )-HT ST  CI 2 1 J ) * ( TEM-A G ( 1 1 > > ) / C AG C II ) - 
1A6CI2)  ) 
IFCHROUC J>-HT)  30«%0«40 
30  CONTINUE 

40  IFCHROyC J)-HT)  50t6Q»60 
50  SITE=120. 
RETURN 

60  SITE=CSI C Jl)*CHT-HROy CJ2> )-SICJ2)*CHT-HR0UCJl) ) )/CHROW( J1)-HR0UC J 
1)  ) 


RETURN 
END 

TAJAPKL  5  PAGES  PRINTED, 

?3IJ08.SEPARAT0R«<? 
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Appendix  III 


Example  Run 

A.  Data 

B.  Program  Execution 

C.  Output 


^DATA.L  spruceoata. 

Data  9R1  sltatq  08/15/80  0*5:23:57  (i) 

1.  1037070016121386     WILLIAMS  CREEK 

2.  1  037070016122020   QA?067'?7   00530300     01032   OIX  6A0111 

3.  103707001612AOP055031 1 1 1010000000002  01 

4.  1 037P700l6l26F16COOflGl37742300203621025S068W06  115 


5. 

1037070016128  0016 

6. 

103707001612Q00101 

01 093077025 

18039 

008200 

1 

7. 

1 03707001612900102 

010931 28041 

09156 

018255 

1 

8. 

1 0370700161 2900103 

01 093070028 

009200 

1 

9. 

1 0370700161?90C1 

0109304001501 

8 

10. 

1 03707001612900104 

01 093000000 

008  029 

0  00 

2 

n  • 

1C3707001612900105 

01093000001 

0  00 

2 

12. 

103707001612900106 

01093000004 

0  00 

2 

13. 

1 03707001612900201 

01093159058 

08168 

019227 

1 

14. 

1 03707001612900202 

01093182057 

06224 

019227 

1 

15. 

1 03707001612900203 

01093058028 

05090 

008355 

2 

16. 

103707001612900204 

01093144055 

06220 

009227 

1 

17. 

103707001612900205 

01093140045 

029227 

1 

18. 

103707001612Q00206 

01 093080033 

019300 

2 

19. 

1 03707001612900207 

01093027013 

02804051 

09300 

2 

20. 

] 03707001612900208 

01093040019 

02809065 

09300 

2 

21. 

1 0370700161 2900209 

01093033016 

09300 

2 

22. 

103707001^12900210 

01 093000003 

0  00 

2 

23. 

1 03707001612900301 

06093161024 

3 

24. 

1 0370700161290C302 

011021 18025 

8 

2 

25. 

103707001612900303 

01102178027 

R 

2 

26. 

1 03707001612900304 

061021 16021 

3 

27. 

103707001612900305 

01093071027 

009300 

1 

28. 

1 03707001612900306 

01093056024 

008300 

2 

29. 

103707001612900307 

01  093197046 

05189 

248252 

1 

30. 

103707001612^00308 

01093186048 

148252 

1 

31. 

1 03707001612900309 

01 0^3065028 

019255 

1 

32. 

1  r.37070016l?9003l0 

01093025012 

08472 

2 

33. 

1 0370700161290031 1 

01093028012 

07451 

2 

34. 

103707001612900312 

01093026009 

07451 

2 

35. 

1  0370700161 ?90C313 

01093015008 

07451 

2 

36. 

1  n37070016l2Q00314 

01 09301 0007 

07451 

2 

37. 

1C3707001612900401 

01093090024 

487228 

2 

38. 

103707001612QC0402 

01093131028 

029255 

1 

39. 

1 03707001612900402 

990931 31 028 

068 

9255 

40. 

1 0370700161 2900501 

010931 17034 

008200 

2 

41  . 

1  03707001 ^12900502 

01 093071030 

017255 

2 

4?. 

103707001612900503 

01093228075 

06200 

007255 

1 

43. 

1 03707001612900504 

01093094033 

008200 

2 

44. 

1 03707001 612900505 

01093013007 

08400 

2 

45. 

1  0370700161 ?900506 

01093000001 

0  00 

2 

46. 

103707001612900507 

01093000000 

0  00 

2 

47. 

103707001 612900508 

01093000000 

0  00 

2 

4ft. 

1 03707001612900501 

990931 17034 

048 

8200 

6 


49.  U}3707001f.l?900503  990932P807S  POO  7255  ? 

END  DATA,   errors:   NONE.   TI«E:      P.7P3  SEC.    IMARE  COUNT:  49 


PRUNtJ/S   R2TMRR, 1 10260241 1      ♦ I N V AOP ♦ 1 » 75/2 0 0 0 

PASG,A  SPRUCEDATA. 

PASG.A  SPRUCEBEETLE. 

9USE     7»  SPRUCEDATA. 

PXOT  SPRUCEBEETLE. XOT 

SPRUCE  BEETLE  STAND  RATING  FOR  386  WILLIAMS  CREEK 
ALEXANDERS   SITE    INDEX   r  50.9557 

TOTAL   TREES  PROCESSED   IN  STAND  =  39 
TREES  PROCESSED   WITH   PBH   >2.6   IN,    =  28 
TREES   PROCESSED   WITH  nRH>10.0   IN.    =  12 

STAND  RISK   RATING   FOR   SPRUCE   BEETLE   IS  MEDIUM 

PFIN 
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Agriculture— CSU,  Fort  Collins 


Rocky 
Mountains 


Great 
Plains 


U.S.  Department  of  Agriculture 
Forest  Service 

Rocky  Mountain  Forest  and 
Range  Experiment  Station 

The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 

RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 

Albuquerque,  New  Mexico 

Bottineau,  North  Dakota 

Flagstaff,  Arizona 

Fort  Collins,  Colorado* 

Laramie,  Wyoming 

Lincoln,  Nebraska 

Lubbock,  Texas 

Rapid  City,  South  Dakota 

Tempe,  Arizona 


•Station  Headquarters:  240  W.  Prospect  St.,  Fort  Collins,  CO  80526 


